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Transverse momentum distributions in ultra-relativistic heavy ion collisions carry considerable information 
about the dynamics of the hot system produced. Direct comparison with the same spectra from p + p colli- 
sions has proved invaluable to identify novel features associated with the larger system, in particular, the "jet 
quenching" at high momentum and apparently much stronger collective flow dominating the spectral shape at 
low momentum. We point out possible hazards of ignoring conservation laws in the comparison of high- and 
low-multiplicity final states. We argue that the effects of energy and momentum conservation actually dominate 
many of the observed systematics, and that p + p collisions may be much more similar to heavy ion collisions 
than generally thought. 



I. INTRODUCTION 



A. Heavy Ion Physics: Relying on Comparison 



The physics program at the Relativistic Heavy Ion Col- 
Uder (RHIC) at Brookhaven National Laboratory is remark- 
ably rich, thanks to the machine's unique ability to collide 
nuclei from 'H to '^^Au, in fully symmetric (e.g. Au+Au or 
p + p ) to strongly asymmetric (e.g. d+An) entrance chan- 
nels, over an energy range spanning more than an order of 
magnitude. The capability to collide polarized protons pro- 
vides access to an entirely new set of fundamental physics, 
not discussed further here. 

Achieving the primary aim of RHIC- the creation and 
characterization of a color-deconfined state of matter and its 
transition back to the confined (hadronic) state- requires the 
full capabilities of RHIC. In particular, comparisons of par- 
ticle distributions at high transverse momentum (pr) from 
Au +Au and p + p collisions, probe the color-opaque nature 
of the hot system formed in the collisions 111 B |3|] . Compari- 
son with reference d+A collisions were necessary to identify 
the role of initial-state effects in the spectra ^. Comparing 
anisotropic collective motion from non-central collisions of 
different-mass initial states (e.g. Au +Au versus Cu + Cu) [J5J] 
tests the validity of transport calculations crucial to claims of 
the creation of a "perfect liquid" at RHIC la]. Indeed, a main 
component of the future heavy ion program at RHIC involves 
a detailed energy scan, designed to identify a predicted critical 
point in the Equation of State of QCD ||711. 

The need for such systematic comparisons is not unique 
to RHIC, but has been a generic feature of all heavy ion 
programs jsl |91], from low-energy facilities like the NSCL 
(Michigan State), to progressively higher-energy facilities 
at SIS (GSI), the Bevatron/Bevalac (Berkeley Lab), AGS 
(Brookhaven), and SPS (CERN). The nature of heavy ion 
physics is such that little is learned through study of a single 
system. 
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B. Bigger is better 

Despite the necessary attention to smaller colliding part- 
ners, these comparisons are ultimately aimed at identifying 
novel aspects of collisions between the heaviest ions, in which 
a highly excited bulk system might be created, with a sufficient 
number of degrees of freedom such that it may be described 
thermodynamically- e.g. in terms of pressure, temperature, 
energy density, and an Equation of State (EoS). If the energy 
density of this system is sufficiently large (typically estimated 
at eciit ~ 1 GeV/fm^ |6]) and its spatial extent considerably 
larger than the color-confinement length ~ 1 frri, then a new 
state of matter- the quark-gluon plasma (QGP) JlOll - may be 
created. Microscopically, such a state might be character- 
ized by colored objects (or something more complicated 11 111 ); 
macroscopically, it represents a region on the phase diagram 
in which the EoS is distinctly different than for the hadronic 
phase lfl2ll . 

Ultra-relativistic collisions between the heaviest nuclei en- 
joy the additional advantage that finite-size effects are small, 
due to high-multiplicity final states. In a small system (e.g. 
final state of an p + p collision) a statistical analysis of 
yields requires a canonical treatment, due to the conserva- 
tion of discreet quantum numbers such as baryon number and 
strangeness Ill3ll . For larger systems, a grand canonical treat- 
ment is more common [e.g.|l4[]. with finite quantum-number 
effects absorbed into, e.g. "saturation factors" (ISJ. 

Due to the large available energy -/s and final-state multi- 
plicity, energy and momentum conservation effects on kine- 
matic observables (spectra, momentum correlations, elliptic 
flow) are generally small. They are accounted for with correc- 
tion factors llm. \V7\\ or neglected altogether. 



C. Multiplicity evolution of single-particle spectra 



Detailed single-particle spectra (e.g. d^N/dpj) have been 
measured at RHIC, for a variety of particle types. Often, the 
shape of the "soft" {pj < 2 GeV/c) part of the spectrum is 
compared to hydrodynamic calculation s lUal or fitted to sim- 
ple "blast-wave" parameterizations [e.g.[l3l to extract the col- 
lective flow of the system. The "hard" sector {pr > 4 GeV/c) 
is assumed to be dominated by the physics of the initial-state. 



high-2^ parton collisions and resulting jets. The physics of 
the "firm" sector (2 < pr < 4 GeV/c) may be the richest of all, 
reflecting the dynamics of the confinement process itself 112011 . 

We would like to focus not so much on the single-particle 
spectra themselves, but on their multiplicity dependence. 
Much has been inferred from this dependence. In the soft 
sector, blast-wave fits to spectra from high-multiplicity final 
states (associated with central A +A collisions) indicate strong 
collective radial flow; the same fits to low-multiplicity final 
states- including minimum bias p + p collisions- appear to 
indicate much weaker flow luHl . This seems to confirm a 
common assumption that p + p collisions are not sufficiently 
"large" to develop bulk collective behaviour 

In the hard sector, one of the earliest and most exciting ob- 
servations IJ, 1221 at RHIC was that the high-pr yield from 
high-multiplicity Am + Am collisions was suppressed, relative 
to appropriately scaled lower-multiplicity A +A or minimum 
bias p+ p collisions. This has been taken as evidence of en- 
ergy loss of hard-scattered partons through a very color-dense 
medium. Meanwhile, the high-pj- part of the spectrum from 
high-multiplicity p + p collisions appear enhanced relative to 
low-multiplicity p + p collisions ll23ll . again suggesting that a 
color-dense bulk system is not produced in p +/9 collisions. 

In this paper, we discuss the effects of energy and 
momentum conservation on the multiplicity evolution of 
single-particle spectra at RHIC. Energy and momentum 
conservation-induced constraints (EMCICs) 16411 have been 
largely ignored in the analyses just mentioned, probably due 
to two reasons. The first is the field's usual focus on the 
highest-multiplicity collisions, where such effects are as- 
sumed small; it seems natural to compare analyses of such 
systems to "identical" ones of smaller systems, forgetting that 
EMCIC effects play an ever-increasing role in the latter case. 
Perhaps the more important reason is that EMCICs do not 
generate "red flag" structures on single-particle spectra; this 
is in contrast to multi-particle correlation analyses, in which 
conservation law-induced correlations may be manifestly ob- 
vious and have even been used to estimate the number of un- 
measured neutral particles in high energy collisions Il24l1 . Es- 
pecially with the enhanced attention on precision and detail 
at the SPS and RHIC, there has been increasing discussion 
of EMCIC effects in 2-particle HIH, 3-particle iHl, and 
A^-particle 112711 observables. Below, we show that EMCIC ef- 
fects on single-particle spectra are also significant, and may 
even dominate their multiplicity evolution. 



D. Organization of this paper 

Several authors [e.g. |281] have discussed finite-number ef- 
fects in statistical models, and many numerical simulations 
of subatomic collisions conserve energy and momentum auto- 
matically [e.g.|29l|30]. However, as pointed out by Knoll IBlll . 
our question- to what extent do EMCICs alone explain the 
multiplicity evolution of spectra?- cannot be addressed from 
these simulations themselves, since dynamic and kinematic 
evolution are interwoven in these models. Thus, in SectionHD 
we discuss a formalism based on Hagedorn's generalization 



of Fermi's Golden Rule, in which dynamics and kinematics 
(phasespace) factorize. This leads to a formula for finite- 
number effects on single-particle spectra, due solely to kine- 
matics, for a fixed dynamical ("parent") distribution. 

In Section |III1 we test the extreme ansatz that all of the 
experimentally-measured multiplicity dependence of single- 
particle spectra is due to EMCICs. We will find surprising 
agreement with this ansatz in the soft sector {pj < 1 GeV/c). 
We will discuss that our formalism is on less firm footing, 
conceptually and mathematically, at much higher pj- Nev- 
ertheless, we explore this regime as well. We find that, in 
the hard sector, the data from heavy ion collisions is clearly 
not dominated by EMCICs, though we point out that ignoring 
EMCICs, especially for p+ p collisions, may be dangerous 
even at highpj-. 

In Sections HVl and IVl we summarize and give an outlook 
for future studies. 



II. EFFECTS OF ENERGY AND MOMENTUM 
CONSERVATION ON SINGLE-PARTICLE SPECTRA 

A. A restricted phase space factor 

Changing the size (central versus peripheral ion collisions, 
e + e collisions, etc) and energy of a collision system will lead 
to different measured single-particle distributions, reflecting 
(1) possibly different physical processes driving the system 
and (2) effects due to phase space restrictions. To focus 
on changes caused by the latter, we consider some Lorentz- 

invariant "parent" distribution / [p) = lE^^, driven by some 
unspecified physical process, but M«affected by energy and 
momentum conservation. For simplicity, we assume that all 
particles obey the same parent distribution. 

In the absence of other correlations, the measured sin gle- 
particle distribution is related to the parent according to llm 

~fc{pi)=f{px)^ (1) 

/ {j\%2 d'p.i5 {p'j - m^) f (Pj)) 5^ {Llr Pi - P) 
J [UU d'Pj?> {pj - mj) f ipj)) 54 (l£, ;,, - P) ' 

where A^ is the event multiplicity. The integral in the numer- 
ator of Equation[T]represents the number of configurations in 
which the A^ — 1 other particles counter-balance pi so as to 
conserve the total energy-momentum P of the event, and the 
denominator, integrating over all A^ particles, is a normaliza- 
tion. 

For N > 10 1I25II . one may use the central limit theorem to 
rewrite the factor in Equation[T]as lfl6l [TtJ [25i I27I1 
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where 



(K) 



dpf{p)-pl 



(3) 



are average quantities and we have set the average three- 
momentum {p(p=i,2.3)) = PiJ=i.2.3/N = 0. We stress that what 
appears in Equation[3]is the parent distribution /, not the mea- 
sured one /( . Hence, for finite multiplicity A^, the averages 
{p") are not the measured ones, which we define as 



{p'X^ dpfA,p)-p; 



(4) 



See also the discussion in Appendix IBJ 

Since pr distributions are commonly reported, we would 
like to estimate EMCIC distortions to pj distributions, inte- 
grated over azimuth and a finite rapidity bin centered at midra- 
pidity. As discussed in Appendix lAl for the approximately 
boost-invariant distributions at RHIC lulll . the measured and 
parent pr distributions are related by 



fc (pt) == f (pt) 



exp 




The notation X indicates the average of a X over the rapidity 
interval used; see Appendix lAlfor details. These averages de- 
pend, of course, on pr and should not be confused with global 
averages {X) (Equation |3]l which characterize the parent dis- 
tribution. 

We would also like to emphasize the fact that since Equa- 
tion |5] depends on the energy of the particle (not just momen- 
tum) it becomes clear that the EMCIC effects are larger on 
heavier particles at the same pj- Thus we should expect that 
the proton spectra will be more suppressed than pion spectra. 

In what follows, we find that ignoring the pi/ {pi) term 
does not affect our results, since the numerator is small for 
the narrow rapidity windows used here, and the denominator 
is large. In discussions below, we set this term to zero. 



volume required to describe the data was considered unre- 
alistically large Ool . Using the mean value theorem, Hage- 
dorn II33I1 generalized the theory so that the "physics term" is 
the interaction matrix element, suitably averaged over all final 
states. 

We wish to make no assumptions about the underlying 
physics (represented by /) driving the observed spectrum fc. 
Rather, we wish to quantify the effect of changing the multi- 
plicity A^, which appears in the phase space term. 

In particular, in the following Section, we compare mea- 
sured single-particle spectra for different event classes. 

We postulate that the parent distributions for, say classes 1 
and 2, are the same (f\ ~ f2). By Equation [3] this implies 
{Pij)i ~ {P/j}! = {P/j)- In this case, the only reason that the 
observed spectra differ (fd 7^ fc.2) is the difference in "mul- 
tiplicity" A^i ^ N2', see Section HTCl for a discussion of A^i . 

To eliminate the (unknown) parent distribution itself, we 
will study the ratio of observed pr distributions, which, by 
Equation|5]becomes 



fc.i (pt) 

fc,2 (Pt) 



exp 



+ 



Kx 



1 



{N2-l)Ni 



{Ni 



l)N2 

1 
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{E^)-{Ey^ (E^)-{Ey^ {E^)-{Ey 



where the constant K is discussed at the end of Section HTCl 
As mentioned at the end of Section III Al numerically unim- 
portant terms in p^ have been dropped. 

Naturally, our postulate cannot be expected to be entirely 
correct; one may reasonably expect the mix of physical pro- 
cesses in p + p collisions to differ from those in Au +Au col- 
lisions. Nevertheless, it is interesting to find the degree to 
which the change in single-particle spectra may be attributed 
only to finite-multiplicity effects. We will find that the pos- 
tulate works surprisingly well in some regions, and fails in 
others. As we will discuss, both the success and failure raise 
interesting and surprising possibilities. 



C. Testing the postulate - how to treat the parameters 



B. Straw-man postulate of a universal parent distribution 

Equations [l]|5] are reminiscent of Fermi's "Golden 
Rule" Il32l[33ll . in which the probability for making a particu- 
lar observation is given by the product of the squared matrix 
element and a quantity determined by available phase space. 
The first term represented the underlying physical process. In 
his original statistical model [32], Fermi originally assumed 
it to be a constant representing the volume in which emitted 
particles were produced; this is equivalent to setting f{p) con- 
stant in Equation[T] While surprisingly successful in predict- 
ing cross sections and pion spectra [e.g. l34ll35ll . the emission 



By our postulate, the phase space factor affecting a pr dis- 
tribution is driven by four quantities. Three, (pj), (E^) and 
(£■), characterize the parent distribution, while A^ is the num- 
ber of particles in the final state. In general, increasing any one 
parameter decreases the effect of phase space restrictions on 
the observed distributions. But what should we expect these 
values to be? They should characterize the relevant system in 
which a limited quantity of energy and momentum is shared. 
They are not, however, directly measurable, and should only 
approximately scale with measured values, for at least five 
reasons discussed here. 

Firstly, the energy and momentum is shared among mea- 
sured and unmeasured (neutrals, neutrinos, etc.) particles 



alike so that A^ should roughly track the measured event mul- 
tiplicity Nmeas, but need not be identical to it. Secondly, 
emission of resonances smears the connection between A^ and 
Nmeas', eg- the emission of an omega meson which later de- 
cays into "secondary" particles (co -^ kkk) increments A^ by 
unity, rather than three, as far as other particles are concerned. 
This latter consideration also affects the kinematic parame- 
ters (pj), (E^) and (£). While energy and momentum are, 
of course, conserved in resonance decay, the aforementioned 
quantities, themselves, are not. Thus, one need not expect 
perfect correspondence between the appropriate kinematic pa- 
rameters in Equation |6] and the measured ones. 

Thirdly, even restricting consideration to primary particles, 
it is unclear that all of them should be considered in the rel- 
evant ensemble of particles sharing some energy and mo- 
mentum. In particular, for space-time extended systems in 
high-energy collisions, the momentum extent of characteris- 
tic physics processes (e.g. string breaking) and causality in 
an approximately boost-invariant scenario suggest that rapid- 
ity slices of roughly unit extent should be considered sepa- 
rate subsystems luol . Of course, the total available energy 
in any event is shared among all such subsystems; i.e. the 
midrapidity subsystem in one event will not have exactly the 
same available energy as that in another event. However, such 
fluctuations are to be expected in any case- surely individual 
collisions will differ from one another to some extent. Thus, 
we repeat our interpretation of the four parameters A^, (pj), 
(E^) and (E): they characterize the scale, in energy and mo- 
mentum, of the limited available phasespace to an A^-particle 
subsystem. 

Fourthly, Equations [T]|6] are appropriate for fixed A^, while 
we will be comparing to measured spectra selected by mea- 
sured charged-particle multiplicity. Thus, A^ would inevitably 
fluctuate within an event class, even if we could ignore the 
above considerations. Naturally, high multiplicity events con- 
tribute to spectra more than low multiplicity events. Simi- 
larly, the average multiplicity in two-particle correlations is 
even more shifted to higher multiplicities. 

Fifthly, as already mentioned in Section HTAl the kinematic 
parameters (pj), (E^) and (E) correspond to the parent dis- 
tribution, which will only correspond identically to the mea- 
sured one in the limit of infinite multiplicity (i.e. no EMCIC 
distortions). See also the discussion in AppendixlB] 

For all of these reasons, we will treat A^, (pj), (E^) and 
(£) as free parameters when testing our postulate against data. 
Our aim is not to actually measure these quantities by fitting 
the data with Equation|6j this is good, since our fits to the data 
only very roughly constrain our four parameters, as discussed 
in the next Section. Rather, our much less ambitious goal is 
to see whether "reasonable" values of these parameters can 
explain the multiplicity evolution of the spectra. 

To get a feeling for these values, we look at p + p colli- 
sions at y/sNN=200 GeV, simulated by the PYTHIA event gen- 
erator (v6.319) |37]. In the model, we can identify primary 
particles, thus avoiding some of the issues discussed above. 
However, the fact that PYTHIA conserves momentum means 
that we access {p'!,)c as defined by EquationH] not the param- 
eters of the parent distribution. Nevertheless, a scale for our 



^max 


(N) 


(pDc 


(pDc 


{E\- 


{E)c 


1.0 


7.5 


0.58 


0.41 


1.45 


0.98 


2.0 


13.4 


0.59 


2.81 


3.89 


1.57 


3.0 


17.9 


0.59 


12.95 


14.01 


2.65 


4.0 


21.5 


0.59 


82.45 


83.55 


5.13 


5.0 


23.4 


0.59 


262.88 


265.03 


8.29 


CO 


23.6 


0.59 


275.23 


276.4 


8.48 



TABLE I: For a given selection on pseudorapidity |r|| < T|,„a.v, 
the number and kinematic variables for primary particles from a 
PYTHIA simulation of p + p collisions at ^/Inn = 200 GeV are given. 
Units are GeV/c or (GeV/c)~, as appropriate. 100k events were used 
and all decays were switched off in simulations. 



^max 


(N) 


(pDc 


{P")c 


{E^}c 


{E}c 


1.0 


16 


0.20 


0.11 


0.40 


0.44 


2.0 


29 


0.21 


0.76 


1.05 


0.68 


3.0 


39 


0.21 


3.5 


3.8 


1.2 


4.0 


47 


0.21 


24 


25 


2.2 


5.0 


51 


0.22 


88 


89 


3.7 



TABLE II: For a given selection on pseudorapidity |ti| < r\max, the 
number and kinematic variables for final state particles (particle in- 
dex KS=1 in PYTHIA ) from a PYTHIA simulation of p + p collisions 
at ^/sjm = 200 GeV are given. 100k events were generated and de- 
fault PYTHIA parameters were used in simulations. Units are GeV/c 
or (GeV/c)-, as appropriate. 



expectations may be set. Table |T] summarizes the result for 
primary particles satisfying a varying cut on pseudorapidity 
where all particle decays where switched off in PYTHIA sim- 
ulations. The results from simulations when resonance de- 
cays were included in simulations are presented in Table |ll] 
These two tables gives us rough estimates of ranges of the to- 
tal multiplicity and kinematic variables that one may expect. 
The bulk component of single-particle spectra is often esti- 
mated with Maxwell-Boltzmann distributions, with inverse 
slope parameters in the range T ~ 0.15 -^0.35 GeV. Again, 
simply for rough guidance, we list Maxwell-Boltzmann ex- 
pectations for our kinematic parameters in Table Hill assuming 
pion-dominated system. 

Finally, a word about normalization- the quantity K which 
appears in Equation |6] Not only energy and momentum, but 
also discrete quantum numbers like strangeness and baryon 





non-rel. limit 


ultra-rel. 
limit 


ifr = 0.15-f0.35GeK 


(pD 


2mT 


8r2 


0.045^0.98 {GeV/cf 


(E^-) 


^T'- + m'- 


I2r2 


0.10H-1.50GeV2 


(E) 


^T + m 


3r 


0.36H- 1.00 GeV 



TABLE III: The average kinematic variables obtained from the 
Maxwell-Boltzmaim distribution /(p) = 4^ --^ e^^'^ using non- 
relativistic and ultra-relativistic limit. A pion gas is assumed. 



number are conserved event by event, affecting the overall 
yield of a given particle species. For example, the related 
phenomenon of "canonical suppression" affects the ratio of 
yields for strange versus non-strange particles, as multiplicity 
varies Il38ll39ll . Since we restrict our attention to energy and 
momentum conservation and the effect on kinematic quanti- 
ties, we are interested in the shape of the spectra ratio, as a 
function of particle momentum, and include a factor K in our 
Equation|6] which should be of order, but not necessarily iden- 
tical to, unity. We do not discuss it further 



ni. TEST OF THE POSTULATE - COMPARISON TO DATA 

We now explore the degree to which the postulate proposed 
above describes the multiplicity evolution of measured pr 
spectra measured in y^s^N — 200 GeV collisions at RHIC. 
As is frequently done, we will separately discuss the "soft" 
iPT < 1 GeV/c) and "hard" (pr > 3 GeV/c) portions of the 
spectra. This separation is not entirely arbitrary, as spectra in 
these two pj ranges are thought to be dominated by quite dif- 
ferent physics, and the multiplicity evolution in the two sec- 
tors is usually interpreted in terms of distinct physics mes- 
sages. 

In the soft sector, the spectral shapes are often consis- 
tent with hydrodynamic calculations [e. g. [isl 14011 . or fitted 
with blast-wave type models [e.g. [T^, ^M, and show evi- 
dence of strong, explosive flow associated with a collective 
bulk medium. This is especially clear in the mass dependence 
of the spectra; the mr (or pr) spectrum of heavy particles 
like protons are significantly flatter than that for pions, in the 
presence of strong flow. The multiplicity evolution in this 
sector suggests that high-multiplicity collisions (say, central 
Au+Au collisions) show much more collective flow than do 
low-multiplicity (say, p + p) collisions ll2lll . Such an inter- 
pretation initially sensible in a scenario in which flow is built 
up through multiple collisions among emitted particles; the 
concept of a collective bulk medium in a very low-multiplicity 
collision is thus usually considered questionable. 

Particle yields at high pj, on the other hand, are generally 
discussed in the context of fragments from high-g' parton 
scatterings in the initial stage of the collision. As the event 
multiplicity in Au +Au collisions is increased, a suppression 
of high-pT" yields is observed, relative to a properly normal- 
ized minimum-bias spectrum from p+p collisions. This sup- 
pression has been attributed to partonic energy loss in the bulk 
medium 61 miHlil. 

The multiplicity evolution of the spectra in p+ p collisions, 
however, shows quite the reverse. Relative to the soft sector, 
the high-pr yields increase as the multiplicity increases; one 
may also say that the pj spectra become less steep as multi- 
plicity increases 112311 . This seems to reinforce the conclusion 
discussed above in relation to the soft sector, that p + p col- 
lisions do not build up a bulk system capable of quenching 
jets. 

Here, we reconsider these conclusions based on the mul- 
tiplicity evolution of the spectra, in light of the phase space 
restrictions discussed above. 



Event selection 


N 


(p2) [(GeV/c)2] 


(£2) [GeV2] 


(E) [GeV] 


p + p min-bias 


10.3 


0.12 


0.43 


0.61 


Au+AulQ-&0% 


15.2 










Au+Au 60-10% 


18.3 










Au+Au 50-60% 


27.3 










Au+Au 40-50% 


38.7 










Au+Au 30-40% 


67.6 










Au+Au 20-30% 


219 










Au+Au 10-20% 


>300 










Au+Au 5-10% 


>300 










Au+Au 0-5% 


>300 




" 





TABLE IV: Multiplicity and parent-distribution kinematic parame- 
ters which give a reasonable description of the spectrum ratios for 
identified particles in the soft sector. See text for details. Note that 
the multiplicity changes with event class; the parent distribution is 
assumed identical. 



A. Soft sector: identified particles in Au +Au versus p + p 

Figure [T] shows mj- distributions for minimum-bias 
p + p collisions and multiplicity-selected Am + Am collisions, 
all at ^/s^ = 200 GeV, reported by the STAR Collaboration 
at RHIC 1I21I1 . For the highest-multiplicity Am +Au collisions 
(top-most filled datapoints), the spectrum for heavier emit- 
ted particles is less steep than the essentially exponential pion 
spectrum. Circles in Figure |2] show the result of fits with a 
blast-wave model lll9ll . They indicate a kinetic freezeout tem- 
perature of about 100 MeV and average collective flow veloc- 
ity about 0.6c for the most central collisions. For lower multi- 
plicity collisions, the freeze-out temperature appears to grow 
to ^ 130 MeV and the flow velocity decreases to ~ 0.25c. 
The STAR collaboration, using a slightly different implemen- 
tation of a blast-wave model, reported essentially identical 
values mi]. 

Ratios of spectra from minimum-bias p + p collisions to 
those from Am +Au collisions are plotted in Figure [3] For the 
filled points, the denominator is the most central Au +Au col- 
lisions, while the open points represent the ratio when the de- 
nominator is from peripheral (60-70% centrality) Am +Au col- 
lisions. Pions, kaons, and protons are distinguished by differ- 
ent symbol shapes. 

The curves show the function given in Equation |6] for the 
kinematic scales given in Table |IV] Clear from the Table is 
that all curves in Figure [3] are generated with the same kine- 
matic variables {p\), (E^) and {E)\ only the relevant multi- 
plicity changes. 

We do not quote uncertainties on the kinematic or multi- 
plicity parameters, as the fitting space is complex, with large 
correlations between them. Furthermore, it is clear that the 
calculated curves do not perfectly reproduce the measured ra- 
tios. However, it is also clear that "reasonable" values of mul- 
tiplicity and energy-momentum scales go a long way towards 
explaining the multiplicity evolution of the spectra, even keep- 
ing physics ("parent distribution") fixed. Our postulate of Sec- 
tion HlB] seems to contain a good deal of truth. 
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FIG. 1: (Color online) Transverse mass distributions for pions (left), kaons (center) and antiprotons (right) measured by the STAR Collab- 
oration for y/sf^[^=200 GeV collisions |2l|l. The lowest datapoints represent minimum-bias p + p collisions, while the others come from 
Au +Au collisions of increasing multiplicity. Filled datapoints are for the top 5% and 60-70% highest-multiplicity Au +Au collisions, and for 
the p + p collisions. 
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FIG. 2: (Color online) Circles show the temperature (top panel) and 
flow (bottom panel) parameters of a Blast-wave model 111911 fit to the 
STAR spectra of Figure [T] as a function of the event multiplicity. 
Squares represent Blast-wave fit parameters to "EMCIC corrected 
spectra," and shaded region represents these results combined with 
systematic errors, as discussed in the text. 



FIG. 3: (Color online) The ratio of the pj distribution from 
minimum-bias p + p collisions to the distribution from 0-5% (filled 
datapoints) and 60-70% (open datapoints) highest multiplicity 
Au+Au collisions; c.f. Figure[T] The ratio of the kaon spectra from 
p + p and 0-5% Au + Au collisions (solid green squares) has been 
scaled by a factor 1.7 for clarity. Curves represent a calculation of 
this ratio (ratio of EMCIC factors) using Equation|6] 



Another way to view the same resuhs is useful. While 
the curves shown in Figure |3] only approximately describe 
the data shown there, one may approximately "correct" the 
measured mj distributions, to account for EMCICs. This is 



shown in Figure IH where the measured min-bias p + p and 
central and mid-peripheral Au +Au spectra have been copied 
from the full points of Figure [T] and are shown by full points. 
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FIG. 4: (Color online) Transverse mass distributions of pions, kaons and antiprotons for minimum-bias p + p collisions and 60-70% and 
0-5% highest multiplicity Au+Au collisions at ^/sNN = 200 GeV. Filled datapoints are the same as in Figure [T] Open triangles represent 
the P + P spectra divided by the lower curves shown in Figure (3] Open circles are the same spectra as the open triangles, except scaled up 
to compare to the spectra from the Au + Au collisions. Open squares represent the spectra from 60-70% highest multiplicity Au+Au events, 
divided by the ratio of upper and lower curves shown in Figure[3] See text for details. 



The open red triangles represent the min-bias p + p spectra, 
divided by Equation |6] with the parameters from Table |IV] 
This "EMCIC-corrected" spectrum is then scaled up to show 
comparison to the spectra from central Au +Au (open red cir- 
cles); the level of (dis)agreement is identical to that between 
the lower datapoints and curves in Figure |3] 

Spectra from the mid-central Au+Au collisions have 
been likewise "corrected." The open squares in Figure |4] 
may be compared to the open circles; again the level of 
(dis)agreement is equivalent to that between the upper data- 
points and curves in Figure [3] 

Spectra themselves contain more information than two- 
parameter fits to spectra. However, much has been made of 
blast-wave fits to measured pr spectra, which suggest a much 
larger flow in central Au +Au collisions, relative to p+ p col- 
lisions. Thus, it may be instructive to see how EMCICs af- 
fect these parameters. In Figure |5] the pr distributions for 
P + P collisions and the six lowest multiplicity selections on 
Au +Au collisions are shown. Blast wave fits to the measured 
spectra, resulting in the parameters shown by red triangles in 
Figure |2] are shown as curves. On the linear scale of the Fig- 
ure, some deviations between the fit and data, particularly at 
the lowest pr for the light particle, is seen. This has been 
observed previously in Blast-wave fits, and may be due to res- 
onances lll9ll46ll . Nevertheless, the fits to measured data are 
reasonable overall, and for simplicity, we do not exclude these 
bins. 

Also shown in Figure |5] are the "EMCIC corrected" spec- 
tra, as discussed above. As already seen in Figure ID these 
differ from the measured spectra mostly for low multiplicity 
collisions and for the heavier emitted particles. Blast-wave 
fits to these spectra are also shown. Especially for the very 
lowest multiplicity collisions, these fits are less satisfactory 
than those to the measured spectra; the "parent distributions" 



extracted via our approximate EMCIC correction procedure 
follow the Blast-wave shape only approximately. Much of the 
deviation is at pr ~ 0.9 GeV/c for protons from the lowest 
multiplicity collisions (upper-right panels). This is the region 
around which the approximations used in deriving the EM- 
CIC correction should start to break down, as discussed in 
Appendix |B] So, two fits are performed: one including all 
datapoints shown (blue squares in Figure |2|, and the other 
excluding proton spectra points with pr > 0.8 GeV/c. The 
resulting range of Blast-wave parameters is indicated by the 
shaded region in Figure |2l There, statistical errors on the fit 
parameters have been multiplied by ^/^x}Jdx)Y. (ranging from 
^ 2 for spectra from p + p collisions to ~ 1 for those from 
mid-peripheral and central Au+Au collisions) and added to 
both ends of the range. Thus, the shaded region should rep- 
resent a conservative estimate of blast-wave temperature and 
flow strengths to the parent distributions. 

In summary, to the extent that the curves in Figure [3] de- 
scribe the ratios shown there- which they do in sign, magni- 
tude and mass dependence, but only approximately in shape- 
the data is consistent with a common parent distribution for 
spectra from all collisions. The residual deviation seen in Fig- 
ure [3] is observed again in different forms in Figures |4] and |2] 
The upshot is that EMCICs may dominate the multiplicity 
evolution of the spectra in the soft sector at RHIC. Extracting 
physics messages from the changing spectra, while ignoring 
kinematic effects of the same order as the observed changes 
themselves, seems unjustified. 

In particular, STAR iH and others fll have fitted the 
spectra with Blast-wave distributions, which ignore EMCIC 
effects. Based on these fits, they concluded that the difference 
in spectral shapes between high- and low-multiplicity colli- 
sions was due to much lower flow in the latter; c.f. Figure |2] 
Recently, Tang et al. |47] arrived to the same conclusion, us- 
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FIG. 5: (Color online) dN /dpf spectra for pions (left), kaons (center), and protons (right) are plotted on a linear scale, as a function of 
event multiplicity. Top panels show spectra for minimum-bias p + p collisions, and the spectra for the six lowest multiplicity selections of 
Au+Au collisions are shown in the lower panels. Filled symbols are the measured data, while open symbols are the "EMCIC corrected" 
distributions, discussed in the text. (For pions, these distributions overlap almost completely.) Blast-wave fits are indicated by the curves. For 
the "EMCIC corrected" spectra, two fits are performed, to estimate systematic errors. The solid line represents a fit to all datapoints, while the 
fit indicated by the dashed line ignores proton yields above pj = 0.8 GeV/c. 



ing a modified Blast-wave fit based on Tsallis statistics. This 
requires introduction of an extra parameter, q, intended to ac- 
count for system fluctuation effects ll48ll . However, contrary 
to the claims in the Tang paper, the Tsallis distribution - with 
or without q - does not account for energy and momentum 
conservation 114911 ; EMCIC effects would need to be added on 
the top of the Tsallis statistics 114911 . Therefore, conclusions 
about flow in low-multiplicity collisions based on these fits 
are suspect. 



An independent measurement of flow would help clarify 
this issue. Two-particle femtoscopy ("HBT") is a sensitive 
probe of collective motion IpOll and has been measured in 
p + p collisions at RHIC Iplll . Any scenario should be able to 
describe simultaneously both the spectral shapes and the mj- 
dependence of the femtoscopic scales. A study of this topic is 
underway. 



B. Soft sector: unidentified particles in multiplicity-selected 
p + p collisions 

While minimum-bias p + p collisions are the natural "ref- 
erence" when studying Au +Au collisions, the STAR experi- 
ment has also measured pr spectra from multiplicity-selected 
p + p collisions [23y. These are reproduced in Figure |6l in 
which the lowest-multiplicity collisions are shown on the bot- 
tom and the highest at the top. Numerical labels to the right 
of the spectra are included just for ease of reference here. 

The solid curve is a power-law fit to the highest-multiplicity 
spectrum (#10), just for reference. This curve is scaled and 
replotted as dashed lines, to make clear the multiplicity evo- 
lution of the spectra. Concentrating on the soft sector for 
the moment, we perform the same exercise as above, to see 
to what extent this multiplicity evolution can be attributed to 
EMCICs. 

In Figure [T] are shown three ratios of spectra, in which the 
second-highest-multiplicity spectrum (#9) is used as the de- 
nominator, to avoid statistical fluctuations associated with the 
highest multiplicity spectrum. Also shown are curves, using 
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FIG. 6: (Color online) Transverse momentum spectra of unidentified 
negative hadrons from p-\-p collisions at ^/Jnn = 200 GeV by the 
STAR Collaboration 12311 . The lowest (highest) dataset corresponds 
to the lowest (highest) multiplicity collisions. The solid line is in- 
tended only to guide the eye and show the shape of the spectrum 
for the highest multiplicity selection. It is rescaled and redrawn as 
dashed lines below, to emphasize the multiplicity evolution of the 
spectrum shape. 



Equation |6] with the energy-momentum scales given in Ta- 
ble |V] 

The spectra reported by STAR are for unidentified nega- 
tive hadrons. In calculating these curves, we assumed that 
all particles were pions. This matters, since the energy terms 
in Equation |6] require the particle mass. We expect the 
energy-momentum scales listed in Table |V] to be affected by 
this simplistic assumption. Particle-identified spectra from 
multiplicity-selected p + p collisions would be required, to do 
better. Given this, and the only semi-quantitative agreement 
between the calculations and measured ratios shown in Fig- 
ure |2l we conclude only that the EMCIC contribution to the 
multiplicity evolution of low-z^j- spectra in p + p collisions is 
at least of the same order as the observed effect itself. 
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TABLE V: Multiplicity and parent-distribution kinematic parame- 
ters which give a reasonable description of the spectrum ratios for 
unidentified particles in the soft sector from multiplicity-selected 
p + p collisions. See text for details. Note that the multiplicity 
changes with event class; the parent distribution is assumed identical. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
p^ [GeV/c] 

FIG. 7: (Color online) Ratio of the pr spectra shown by full points 
in Figure |6] Spectra for the lowest-multiplicity (red triangles), fifth- 
lowest (green triangles) and seventh-lowest (squares) multiplicity 
collisions are divided by the spectrum for the second-highest mul- 
tiplicity collisions. Curves represent a calculation of this ratio (ratio 
of EMCIC factors) using Equation[6l see text for details. 



C. Segue: From the soft to the hard sector 

Figure[3]shows the central result of this paper: namely, that 
the multiplicity evolution of the mass and pr dependence of 
single particle spectra in the soft sector may be understood 
almost entirely in terms of phase-space restriction with de- 
creasing event multiplicity. 

Plotted in that figure is the ratio of spectra from low- 
multiplicity events over spectra from high-multiplicity events. 
Experimental studies sometimes show this ratio's inverse, of- 
ten called Raa 0211 ■ While of course the same information is 
shown in both representations, we choose that of Figure |3] for 
two reasons. The first is to emphasize the effects of EMCICs, 
the topic of this paper; these are, generically, to suppress the 
particle yield at high energy and momentum, particularly for 
low-N final states. (In multiparticle distributions, they also 
generate measurable correlations 12511 .) 

The second reason is to stress that we have been discussing 
spectra in the soft sector, whereas the ratio Raa is generally 
studied at high pj. At large pj-, we expect that a purely 
EMCIC -based explanation of the multiplicity evolution of the 
spectra might break down, for two reasons. Firstly, even if 
particles of all momenta shared phase-space statistically, our 
approximation of Equation |2] is expected to break down for 
energies much above the average energy, as discussed in Ap- 
pendix |B] Secondly, it is believed that the high-/77 yield has 
a large pre-equilibrium component; thus, high-pj- particles 
might participate less in the statistical sharing of phase-space, 
as discussed in Section HTCl 

As we discuss in the next Section, EMCICs surely do 
not dominate the multiplicity evolution of the hard sector in 
heavy ion collisions. For interpreting high-pj- spectra from 
multiplicity-selected p + p collisions, accounting for EMCICs 
may or may not be important. In order to make the connection 
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FIG. 8: (Color online) The same data and curves as shown in Fig 
ure|7l but plotted over the entire measured pj range. 



to Figure [3] we will plot spectra from low-multiplicity colli- 
sions over those from high-multiplicity, as well as the inverse, 
to make the connection to Raa- 



D. Spectra In the hard sector 

The generic effect of EMClCs is to suppress particle yields 
at energy-momentum far from the average value. The effect 
is stronger for lower multiplicity A^. It is clear, then, that EM- 
ClCs cannot account for the multiplicity evolution of the spec- 
tra at high pt in Au+Au collisions, since high-multiplicity 
collisions are observed to have more suppression at high pr 
than do low-multiplicity collisions yll. Thus, we conclude 
that our postulate fails for Au +Au collisions at high pr; the 
"parent distribution" describing the underlying physics in this 
region does, indeed, change with multiplicity. 

But in p + p collisions, the multiplicity evolution in the 
hard sector is opposite to that in Au +Au collisions. In partic- 
ular, in p+ p collisions, the yield at high pr (relative to lower 
Pt) is increased as multiplicity increases, as is clear from Fig- 
ure |6j similar results have been observed m'p + p collisions 
at the Tevatron |52], ISR [53], and SppS H. A "hardening" 
of the spectrum with increasing multiplicity goes in the same 
direction as would EMCIC effects. To what extent can EM- 
ClCs account for the multiplicity evolution of spectra from 
p + p collisions, in the hard sector? 

Some insight on this question may be gained from Figure[8] 
in which the data and curves shown in Figure|2]are plotted out 
to Pt ^6 GeV/c. Clearly, the calculated suppression function 
(Equation|6l) fails dramatically at high pT- 

We recall that Equations |2] and |6] are based on the cen- 
tral limit theorem (CLT), which naturally leads to Gaussian 
distributions. As discussed in Appendix iBl one expects the 
breakdown of the CLT approximation in the far tails of the 
distribution- e.g. when pj ^ (Pt)- Thus, any inferences we 
make about EMCIC effects in the hard sector remain quali- 
tative. Nevertheless, the level of disagreement between the 



FIG. 9: (Color online) "Rpp," the analogue of "i^cf " used in heavy 
ion collisions. The spectrum from the highest-multiplicity p + p col- 
lisions are divided by spectra from lower-multiplicity collisions (see 
filled datapoints in Figure [S^. The data and curves are simply the 
inverse of those shown in Figure [8] 



calculations and measurements leads us to conclude that EM- 
ClCs do not fully explain the multiplicity evolution of pT 
spectra in p + p collisions in the hard sector. 

However, this, in itself, raises a fascinating possibility. Fig- 
ure[8]shows that, relative to high-multiplicity p + p collisions, 
the suppression of high-pT yields from low-multiplicity col- 
lisions is not as strong as one expects from our simple postu- 
late. Said another way, the high-pT "enhancement" in high- 
multiplicity collisions may not be as large as one expects from 
phasespace considerations alone. This is emphasized in Fig- 
ure|9] in which is plotted "Rpp\ the ratio of the spectrum from 
high-multiplicity to lower-multiplicity collisions; Rpp is the 
analog of Rcp from heavy ion collisions [Si] . 

The motivation for studying quantities like Raa and Rqp 
(and now Rpp) is to identify important differences between 
one class of collisions and another. Presumably, one is inter- 
ested in physics effects (jet quenching, etc.), above and be- 
yond "trivial" energy and momentum conservation. Thus, it 
makes sense to attempt to "correct" for EMCICs by divid- 
ing them out as we did in Section ITlIAl keeping in mind the 
caveats just discussed. 

The result of this exercise is shown in Figure (TO] in which 
the datapoints from Figure |9] are divided by the curves from 
the same Figure, to form a new quantity, R' . Explicitly, the 
green circles on Figure[TOl which compare multiplicity selec- 
tions #9 and #4 are given by 




where the relevant quantities from Table [V] are used. Again, 
all particles are assumed to have pion mass. Qualitative 
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FIG. 10: (Color online) Rpp (c.f. FigurelU divided by the EMCIC 
contribution to Rpp, as calculated by Equation [7] 



though it is, Figure [TOlraises the possibility that, when "triv- 
ial" EMCICs are accounted for, the high'/jj- yield from high- 
multiplicity p+ p collisions is suppressed relative to low- 
multiplicity collisions, a trend in the same direction as that 
observed in Au +Au collisions. 

In the hard sector, our estimates are mathematically and 
conceptually too simplistic to decide whether this implies "jet 
quenching" in high-multiplicity p + p collisions. However, 
it is quite clear that conservation-induced phasespace restric- 
tions might be sufficiently large in the hard sector, so that a 
high-pT" "enhancement" in high-multiplicity p + p collisions 
turns into a "suppression," when these effects are accounted 
for. Extracting physics messages (e.g. about mini-jet pro- 
duction or jet quenching) from the multiplicity evolution of 
p+p spectra is a non-trivial task, in light of this potentially 
huge background effect. At the very least, EMCICs should 
not be ignored, as they usually are, when extracting physics 
messages. 



IV. SUMMARY AND DISCUSSION 

The study of relativistic heavy ion collisions is, by its very 
nature, heavily dependent on comparative systematics. Phys- 
ical models or hypotheses are most stringently tested when 
predictions for a given observable are compared to measure- 
ments for a range of global collision conditions. Even aside 
from specific models, much qualitative information may be 
gleaned simply through study of the evolution of an observ- 
able as collision conditions- quantified by global variables- 
change. 

Since the goal is to probe an interaction or transition charac- 
terized by a dimensionful scale (confinement length ~ 1 fm), 
perhaps the most important global variable is event multiplic- 
ity, which on average reflects the size of the system generated 
in the collision. 

Directly measurable is the multiplicity evolution of exper- 
imental observables. This evolution is driven by (1) the evo- 



lution of the underlying physics- which is of direct interest 
and (2) kinematic phase-space restrictions (EMCICs)- which 
are presumably less interesting. It may be hazardous to ignore 
the latter effect and make inferences on the former, particu- 
larly since phase-space restrictions have an obvious explicit 
multiplicity dependence. In this study, we have quantitatively 
estimated the degree to which phase-space restrictions may 
affect physics inferences based on measured data. 

We have focused on the multiplicity evolution of single par- 
ticle spectra. In previous published studies, analyses which 
have ignored EMCICs have inferred much from this evolu- 
tion. In particular, there have been conclusions that spectra 
from central Au+Au collisions exhibit greater collective ra- 
dial flow than do those from peripheral Au +Au or p + p colli- 
sions. Using an expression to approximately account for EM- 
CIC effects, we have shown that the multiplicity evolution of 
the spectra may be dominated by such effects, rather than any 
change in the underlying physics. 

In particular, we have tested the extreme postulate that the 
driving physics, characterized by a parent distribution, is iden- 
tical for p + p collisions and Au +Au collisions of all central- 
ities. Since the parameters characterizing the parent distribu- 
tion and the system multiplicity A^ were fitted, our test is not 
perfect. Some multiplicity evolution of the parent distribution 
itself may exist, and may not be easily separable from EM- 
CICs. Our point is that, with "reasonable" parameters, much 
of the data systematics is readily understood in terms of a uni- 
versal parent distribution in the soft sector, and similar high- 
er yield suppression in p + p and Au +Au collisions. 

In the soft sector [pr < 1 GeV/c) this postulate worked sur- 
prisingly well. The changes in mj distributions, as the colli- 
sion multiplicity is changed, are almost entirely due to EM- 
CICs. "Correcting" the spectra for EMCICs, an approximate 
procedure along the lines of Fermi's Golden Rule, reveals al- 
most universal parent distributions. 

While the spectra themselves carry more information than 
fits to the spectra, it was interesting to find that blast-wave fits 
to the "EMCIC-corrected" spectra show that low multiplicity 
Au+Au collisions, and even p + p collisions, are character- 
ized by very similar flow and temperature values as for spectra 
from Am +Au collisions. This contrasts strongly with previous 
conclusions and assum ptio ns about collectivity in small sys- 
tems. Blast-wave 11191 12111 or modified Blast-wave ll47ll fits 
which ignore EMCICs, may yield unreliable results for low- 
multiplicity final states. 

The same analysis of pj spectra of unidentified hadrons 
from multiplicity-selected p + p collisions yielded similar re- 
sults, though the multiplicity evolution of the spectra was only 
roughly explained by our postulate. This is to be expected, for 
several reasons. Firstly, our approximate expression to ac- 
count for EMCICs was based on the central limit theorem, 
which begins to break down for the very small multiplicities 
involved. Secondly, the lack of particle identification led to 
a simple assumption that all particles were pions. Neverthe- 
less, it was clear that EMCICs can go a long way towards 
explaining the multiplicity evolution of the pr spectra in the 
soft sector 

EMCIC effects on momentum distributions are expected to 
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be large at higher pr, where a single particle may consume 
much of the total available energy. However, the approxima- 
tions behind our EMCIC factor should begin to break down 
at high pT- Unlike our results in the soft sector, we would be 
on shaky ground to draw firm conclusions from our studies in 
the hard sector. Nevertheless, we applied our formalism to ob- 
tain a rough estimate the magnitude of restricted phase-space 
effects at high pj- 

Firstly, we immediately realized that the well-known "high- 
er suppression" for central Au +Au collisions can not be ex- 
plained by EMClCs, as these effects would cause the opposite 
behavior (i.e. "high-pj- enhancement") from what is experi- 
mentally observed. Thus, our postulate fully breaks down at 
high pt- there is a difference in the physics (parent distribu- 
tion) in the hard sector 

Turning to the multiplicity-evolution of pr spectra from 
p + p collisions, however, the measured effect goes in the 
same direction as that expected from EMCIC effects. Still 
keeping in mind the caveats behind our expression at high 
momentum, we estimated that the high-pj- enhancement ex- 
pected from EMCICs should be at least as large as that ob- 
served in the data. Again, we do not conclude, but suggest 
that the multiplicity-evolution of the parent distributions in 
p + p collisions might in fact reveal a high-/77 suppression 
for high multiplicity collisions, reminiscent of the effect mea- 
sured in heavy ion collisions. 



V. CONCLUSIONS AND OUTLOOK 

Our results suggest that the multiplicity evolution of the soft 
portion of the pr spectra in collisions at RHIC is dominated 
by phase-space restrictions. Effects due to actual changes in 
physics (the parent distribution) are subdominant. This sug- 
gests one of two possibilities. 

Firstly, one may take the common assumption that the 
physics underlying the soft particles from A +A and p + p col- 
lisions is quite different, say bulk behavior versus string 
breaking, respectively. In this case, our results suggest that 
single-particle spectra are too insensitive to distinguish very 
different physics scenarios, and physics conclusions (say, ra- 
dial flow in A + A collisions) based on them are questionable. 

On the other hand, the single-particle spectra may well re- 
flect the underlying physics. If energy and momentum con- 
servation effects are taken into account, the low-z^j- spectra 
indicate that p + p collisions display as much collective radial 
flow as do Au +Au collisions. In the larger system, this col- 
lective behavior is usually considered to arise from a (perhaps 
only partially) thermalized bulk system. 

The question naturally arises: isn't it impossible for a sys- 
tem as small as that created ina p + p collision to form even a 
partially thermalized bulk system which develops flow? The 
answer is not obvious. After all, estimates set the timescale 
for complete thermalization in central Au +Au collisions be- 
low 1 fm/c lUSllJOll . via a mechanism that may be driven more 
by fluctuating color fields than by classical rescattering pro- 
cesses 1I55I ,and references therein]. Perhaps the possibility 
that similar processes have sufficient time to thermalize a sys- 



1 fm should not be dismissed out of 



tem on the scale of 
hand. 



Indeed, in the Uterature one finds frequent suggestions 11531 
|5a,|53.|5M|59l|601, based on single-particle spectra, that high 
energy particle collisions generate flowing bulk systems and 
perhaps even Quark-Gluon Plasma; see also the recent review 
by Weiner 16111 . By partially removing the obscuring effects of 
EMCICs, we have more directly compared proton colhsions 
to heavy ion collisions (at the same energy and measured with 
the same detector), for which a flow-based interpretation is 
generally well accepted. 

If a bulk system is created in p+ p collisions, might it 
"quench" jets as the medium does in Au +Au collisions? This 
was, after all, the original proposition of Bjorken 116 211 . The 
signature of such quenching would be a suppression of par- 
ticle yields at high pr in high-multiplicity collisions, relative 
to those at lower multiplicity. While our formalism is insuf- 
ficiently reliable at high pj to draw firm quantitative conclu- 
sions, such a suppression may possibly be present, though ob- 
scured by EMCICs in measured spectra. 

Increased focus on the relationship between large and small 
systems created in ultrarelativistic collisions is called for Ex- 
perimental programs at the Large Hadron Collider will very 
soon open up important avenues in this study. In particular, 
the experiments will measure first p + p collisions at record 
collision energies, with event multiplicities similar to Cu + Cu 
or semi-peripheral Am +Au collisions at RHIC. Soft sector pj 
distributions will likely be among the first observations re- 
ported. Later, with identical acceptance and techniques, the 
same experiments will then measure much larger systems cre- 
ated in Pb + Pb collisions. The direct comparison afforded by 
this data should help answer the question of whether a bulk 
system created in hadronic collisions is qualitatively different 
than that created in collisions between the heaviest ions, or 
merely a smaller version of it. 

The nature of relativistic heavy ion studies depends upon 
comparison of "small" and "large" collision systems, each 
of which may be driven by distinct, non-trivial physics pro- 
cesses. In performing such comparisons, we must not neglect 
the "trivial" effect of energy and momentum conservation, and 
its explicit dependence on collision size. 
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FIG. 1 1 : (Color online) EMCIC factor calculated using the numerical 
averaging of Equation |A2] and the approximation of Eauation lA3l 



1 dN 



FIG. 12: (Color online) j^^ obtained from GENBOD events run for 
the same average energy ({E)c = 1 GeV) but different multiplicities: 
A' = 5, 10, 15, 20, 30, 40 pions. 



APPENDIX A: EMCIC FACTORS FOR RAPIDITY- AND 
ANGLE-INTEGRATED pT DISTRIBUTIONS 

Equation |2] gives the EMCIC correction factor to the triple 
differential spectrum f [p). Experimental measurements often 
report pj distributions integrated over angle and a range of 
rapidity, i.e. 

fciPT) = j- / dif dyfc{p,,Pr,P,,E). (Al) 

In the absence of a triple-differential measurement, we con- 
sider azimuthally-symmetric distributions, and {p^) = {p^) = 
{p\)j2. At midrapidity at RHIC, it is reasonable also to as- 
sume a boost-invariant parent distribution. In this case, only 
part of the EMCIC factor remains in the rapidity integral: 



fcXPT)=f{PT)- 
1 



A^ 



E^ 



N-l 



dyexp 



exp 



9 
-Pf 



{N-l){pl) 



-1 



Pt 



2{N-l)\{pl) 

2E(E) (£)2 



(A2) 



E^)- EV E')- E) 



To arrive at a closed form for our EMCIC factor, we ap- 
proximate the average of the exponential with the exponential 
of the average, i.e. 



fc{PT)^f{PT) 



exp 




pI 



Hn-\)\{pI) (p2) 

£2 2E{E) 



(A3) 



(£2)_(£)2 (£2)_(£)2 (£2)_(£)2 



This expression is reproduced in Equation|5] 
Here, the rapidity-averaged quantities are 



Pi 



1 



2^,, 



pldy^ 



ymax 



£2 = 



2y 



E^dy = mf- 



max J — V 



_ 1 

E 



sinh(23'„ 

^ymax 



/sinh(23VMT) 

\ \^max 



2yn 



/y'liiax 
Edy- 
yiucix 



- nij- 



sinh(3'„ 



(A4) 



(A5) 



(A6) 



The approximation used in going from Equation IA2 1 to 
is well-justified for typical numerical values used in this study. 
Figure [TT] shows a numerical integration of the EMCIC factor 
from Equation IA2I (labeled "exact") and Equation I A3 1 ("ap- 
proximation") for values indicated in the Figure. 



APPENDIX B: REGION OF APPLICABILITY FOR THE 
EMCIC FORMULA 

The exact expression for the phase space integral of Eq. [T] 
was approximated by that in Eq. |2] through an appeal to the 
Central Limit Theorem. Discrepancies between the exact ex- 
pression and the approximate Gaussian functional form will 
become more apparent in the tails of the distribution. For 
example, our approximate phase space suppression function 
never vanishes, thus permitting a tiny but finite probability 
for a particle to carry more energy than that of the entire sys- 
tem! In this Appendix, we perform simple numerical calcu- 
lations with the GENBOD computer program ll63ll . to estimate 
the range of quantitative reliability of Equation |2] 

Given a total energy Eto(, multiplicity A^ and list of particle 
masses, GENBOD produces phasespace-weighted events of A^ 
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FIG. 13: (Color online) Blue points are -^ obtained from GEN- 
BOD events run for A' = 20, {E) = 1 GeV. Black solid curve is an 
exponential, the assumed parent distribution; c.f. Equation IB II Red 
dashed curve is the exponential times the EMCIC factor, as per Equa- 
tionlB2l 



4-momentaby filling Lorentz-invariant phase space according 
to the Fermi distribution. 



f 



d^N 



1 dN 



; 2E — 5" = oc e 

dp^ 2np dE 



-E/i; 



(Bl) 



where ^ characterizes the slope of the energy distributions. 
Since it is (l/p) -dN /dE which is exponential and not (l/p) • 
dN /dE, the inverse slope ^ should not be considered a "tem- 
perature," but only a parameter characterizing the parent dis- 
tribution. 

As a result, generated particles in an event are correlated 
only by energy and momentum conservation. Thus, EMCIC 
effects on the calculated single-particle spectrum, fc (p), are 
given precisely according to Equation[T] 

To evaluate the region of validity of Equation |2] we use 
Eq. IB II as a parent distribution, f (p). Results of this exer- 
cise are presented on Figure [12] which shows energy spectra 
from GENBOD events with the same average energy per par- 
ticle {E)c = Etot/N = 1 GeV, but different multiplicity A^. As 
expected, in the limit of large A^, fc (p) ^ f [p], and it is clear 
that the plotted distribution is increasingly well-described by 
an exponential, as A' increases. 

It is appropriate here to point out why we wish to identify 
the parent distribution in the first place, rather than following 
the procedure outlined in Section HTB] There, the parent dis- 
tribution cancels when taking the ratio of two measured spec- 
tra fc,\/ fc,2, using the postulate that the parent distributions 
/i and /2 are identical. In contrast, the parent distributions 
for the different GENBOD spectra shown in Figure [12] are as- 
suredly yiot the same. Those spectra came from event samples 



having the same {E)c (c.f. Eq.|4]i, and thus different (E) (c.f. 
Eq.[3]l, implying different parents. 

Having at hand a functional form for the GENBOD parent 
distribution, we may test our approximate formula for the 
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FIG. 14: (Color online) Blue points are ^ |^ obtained from GEN- 
BOD events run for A' = 20,{E) = 1 GeV, divided by exp(-£/Q; 
i.e. the blue points from Fig.[T3]divided by the black full curve from 
the same figure. Red dotted line is the EMCIC factor; i.e. the red 
dotted curve from Fig. [13] divided by the black full curve from the 
same figure. 



phasespace modification factor, by fitting the calculated spec- 
trum according to 



dNc 
'dE 



= A ■ p-e 



N 



.,-£/? 



(B2) 



A'-l 



exp 



3p2 



+ - 



2{N-l)J\{p^) 

2E{E) {Ef 



;£2)-(£)2 (£2)-(£) 



+ ■ 



'E^ - (EV- 



where we used the fact that GENBOD generates particles 
isotropically so that < p^ >=< p^ >~< p} >= j < p^ >. 
Since A^ is a known quantity, and {E), (E^) and (p^) may be 
directly calculated from ^, the fit of Equation lB2l has only two 
parameters: the overall normalization A, which is unimportant 
to us, and ^, which characterizes the parent distribution. 

The results are shown in Figure [13] and, for better detail, in 
Figure [14] For the case here, which is typical of that in the 
data, we see that our approximation begins to break down for 
particle energies E >2^3 (E) . Above this range, our approx- 
imation (e.g. Equation|6]) should only be taken qualitatively. 
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